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ABSTRACT 

A  global,  six-level  primitive  equation  model  was  numeri- 
cally integrated  to  examine  the  interaction  of  a  planetary- 
wave  with  a  growing  mid-latitude  baroclinic  wave.   The 
purpose  of  this  study  was  to  determine  the  contribution  of 
the  baroclinic  wave  to  the  initiation  of  a  monsoon  surge. 
Three  experiments  were  run.   Experiment  I  depicted  the 
effect  of  the  planetary  wave  on  the  southward  migration  of 
cold  air.   The  cyclogenesis  of  a  mid-latitude  baroclinic 
wave  was  shown  in  Experiment  II.   Experiment  III  combined 
the  growing  baroclinic  wave  with  the  stable  planetary  wave 
in  an  attempt  to  simulate  a  cold  surge.   Characteristics  of 
a  monsoon  surge  were  present  in  the  numerical  simulation  but 
not  to  the  extent  or  detail  of  observed  surges .   The  results 
indicated  the  important  role  of  planetary  wave  advections 
in  the  formation  of  cold  surges . 
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I.   INTRODUCTION 

The  Asian  winter  monsoon  has  been  an  area  of  very  active 
research  over  the  past  few  years.   Winter  monsoons  are 
dramatic  manifestations  of  land-sea-air  interaction.   These 
phenomena  consist  of  a  mixture  of  all  scales  of  motions,  and 
the  various  interactions  are  still  not  well  understood. 
Although  the  fundamental  structure  of  the  monsoon  has 
frequently  been  observed,  the  dynamical  mechanisms  are 
still  unknown.   Recent  areas  of  investigation  include  the 
role  of  the  winter  monsoon  in  mid-latitude  tropical  inter- 
action, monsoonal  variation  and  the  forcing  mechanisms  for 
small  scale  monsoonal  variations. 

The  monsoon  is  a  three-dimensional,  planetary  scale  wind 
regime  which  exhibits  a  strong  seasonal  dependence.   It 
exhibits  two  distinct  seasonal  low-level  wind  flow  patterns-- 
a  winter  outflow  from  a  cold  continental  anticyclone  and  a 
summer  inflow  into  a  continental  heat  low.   The  winter  anti- 
cyclone and  the  summer  heat  low  result  from  land-sea  tempera- 
ture differences  and  the  varying  position  of  the  sun.   The 
shapes  of  the  continents  and  their  topographies  also  play 
a  role  in  producing  the  considerable  regional  and  temporal 
variability  of  the  monsoons. 

Though  centered  far  to  the  north,  the  Siberian  anti- 
cyclone dominates  the  Southeast  Asian  winter  (Ramage,  1971). 
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The  low  heat  capacity  of  the  land  relative  to  the  open  ocean 
to  the  east  and  south  ensures  the  development  of  the  cold 
anticyclone  with  surface  flow  from  land  to  ocean.   A  com- 
pensating flow  from  ocean  to  land  occurs  in  the  upper 
troposphere.   The  Himalayas  block  the  southward  movement  of 
the  extremely  cold  air  from  the  Siberian  anticyclone.   As  the 
anticyclone  becomes  established  over  Siberia,  a  surface 
cyclone  develops  over  Indonesia  and  Malaysia  as  a  major 
planetary  scale  convective  area  is  shifted  eastward  from  its 
summer  position  near  India  (Baker,  1983). 

Chang  and  Lau  (1980)  analyzed  the  monthly  mean  velocity 
potential  field  at  200  mb  for  December  1974  and  showed  a 
broad  region  of  divergence  over  Indonesia.   This  broad  area 
of  divergence  is  indicative  of  deep  rising  motions  and  large 
scale  overturning  in  the  atmosphere.   A  region  of  maximum 
convergence  over  northern  China,  in  the  vicinity  of  the 
Siberian  anticyclone,  was  also  noted. 

A  strong  Hadley  cell  dominates  the  winter  tropospheric 
circulation  of  this  region.   This  thermally  direct  circu- 
lation involves  the  descent  of  cold  air  over  the  Siberian 
anticyclone,  the  southward  flow  of  air  at  the  surface,  the 
ascent  of  warm  air  over  Indonesia  and  a  northward  return 
flow  at  200  mb. 

Although  the  large  scale  features  of  the  monsoon  are 
identifiable  from  year  to  year,  there  is  large  interannual 
variability  in  strength  and  location.   Velocity,  temperature 
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and  pressure  fields  vary  on  short  time  scales.   One  such 
variation  is  the  monsoon  surge,  also  known  as  the  "cold 
surge."   These  events  occur  at  intervals  of  several  days 
to  two  weeks  during  the  winter  monsoon  (Lim  and  Chang,  1981) . 

Observational  studies  by  Danielson  and  Ho  (1969),  Chang 
et  al.  (1979),  Murakami  (1979),  Chang  and  Lau  (1980,1982) 
and  Chang  et  al.  (1983)  discuss  the  characteristics  and 
effects  of  the  monsoon  surge  at  the  synoptic  and  planetary 
scales.   The  characteristics  and  effects  described  in  these 
studies  are  summarized  in  the  following  paragraphs. 

Over  northern  Asia  there  is  intense,  deep  cold  air 
advection  associated  with  the  passage  of  a  mid-tropospheric 
trough.   The  trough  intensifies  rapidly  over  northern  Japan. 
The  baroclinicity  increases  with  the  anticyclogenesis  over 
Siberia  and  cyclogenesis  in  the  strong  baroclinic  zone  off 
the  east  China  coast.   As  the  anticyclone  intensifies, 
subsidence  increases  which  strengthens  the  east  Asian  Hadley 
circulation.   This  intensification  increases  the  upper-level 
ageostrophic  flow,  causing  an  acceleration  of  the  east  Asian 
subtropical  jet  near  Japan. 

The  initiation  of  the  monsoon  surge  occurs  when  the  upper- 
level  trough  moves  eastward  and  the  pressure  gradient  along 
the  east  China  coast  increases.   When  the  pressure  gradient 
force  exceeds  the  Coriolis  force,  the  normal  geostrophic 
balance  is  broken,  resulting  in  cross-isobaric  flow  from 
high  to  low  pressure.   This  accelerating  northeasterly 
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low-level  flow  toward  lower  pressure  is  the  monsoon  surge. 
The  surge  reaches  the  equatorial  South  China  Sea  in  12-24 
h.   In  the  equatorial  regions,  convection  intensifies  within 
one  to  two  days,  sustaining  the  already  enhanced  Hadley 
circulation.   The  upper-level  outflow  from  the  convection 
accelerates  eastward  and  westward  along  the  equator  as  well 
as  northward.   The  surge  ends  when  the  mid-latitude  trough/ 
ridge  pattern  moves  far  enough  east,  thus  diminishing  the 
strong  subsidence  over  the  anticyclone. 

The  monsoon  surge  has  two  stages  separated  by  a  few 
hours  to  one  day  depending  on  the  location  of  the  observing 
site  (Lim  and  Chang,  1981;  Chang  et  al. ,  1983) .   The  pres- 
surge  surge  is  the  first  stage.   It  is  the  leading  edge  of 
the  downward  accelerating  air  with  a  propagation  speed 
characteristic  of  internal  gravity  waves.   A  frontal  passage 
is  the  second  stage  which  moves  with  advective  speeds  and 
is  defined  by  a  sharp  decrease  in  surface  dewpoint. 

The  cause  or  nature  of  monsoon  surge  forcing  has  been 
examined  by  Lim  and  Chang  (1981)  and  Baker  (1983).   Linearized 
shallow-water  equations  on  an  equatorial  beta-plane  were 
used  by  Lim  and  Chang  (1981)  in  an  analytical  study  of 
monsoon  response  to  simulated  anticyclogenesis .   They  sug- 
gested that  the  large  features  of  synoptic  scale  motions 
in  the  east  Asian  winter  monsoon  could  be  explained  by 
pressure  forcing  using  equatorial  beta-plane  dynamics.   No 
consideration  was  given  to  orography,  a  variable  mean  zonal 
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wind,  the  vertical  structure  of  monsoon  surges  or  the  upper 
tropospheric  response  to  the  surge. 

Baker  (1983)  used  a  global,  six-level  primitive  equation 
model  to  examine  the  interaction  of  a  mid-latitude  baro- 
clinic  wave  with  topography.   A  well-developed  baroclinic 
wave  could  initiate  a  monsoon  surge,  however,  the  surges 
were  weak  with  limited  extent.   Baker  (19  83)  concluded  that 
other  forcing  mechanisms  were  required  to  simulate  the 
surges  observed  during  the  winter  monsoon. 

Shaffer  et  a_l.  (1984)  have  shown  that  cold  surge  events 
are  frequently  preceded  by  rapid  growth  in  the  planetary 
waves.   It  appears  that  the  amplified  planetary  wave  helps 
advect  the  cold  air  toward  the  equator.   The  purpose  of  this 
study  is  to  examine  the  interaction  of  a  planetary  wave 
with  a  growing  baroclinic  wave.   The  simulations  will  be 
carried  out  with  a  global,  six-level  primitive  equation 
model.   The  same  baroclinic  jet  will  be  used,  as  was  em- 
ployed by  Baker  (1983).   The  baroclinic  wave,  which  is  unstable 
with  the  jet,  will  have  wave  number  8.   The  planetary  wave 
will  be  represented  by  wave  number  4 ,  which  is  stable  with 
the  jet. 
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II.   DESCRIPTION  OF  MODEL 

A  version  of  the  UCLA  potential  enstrophy  conserving 
general  circulation  model  is  used  in  these  experiments.   The 
atmosphere  is  adiabatic,  dry  and  f rictionless .   The  governing 
equations  of  the  model  are  the  primitive  equations  (listed 
in  Appendix  A) .   The  model  design  and  dynamic  processes 
are  described  by  Arakawa  and  Lamb  (1977,1981)  .   The  version 
used  for  these  experiments  is  summarized  below. 

A.   THE  HORIZONTAL  GRID 

Arakawa ' s  scheme  C  is  used.   The  prognostic  variables- 
wind  (u,v),  potential  temperature  (a)  and  terrain  pressure 
.  (77) — are  staggered  horizontally  as  shown  in  Figure  2.1.   This 


Tr.e 


Figure  2.1. 


Horizontal  distribution  of  variables  in 
the  UCLA  model  grid  (Arakawa  and  Lamb, 
1977)  . 


grid  arrangement  requires  more  computer  time  than  an  unstag- 
gered  grid,  but  it  provides  a  much  better  structure  for  the 
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shorter  waves  arising  from  geostrophic  adjustment,  smaller 
scale  heating,  nonlinear  interaction  or  topographic  effects. 
Grid  resolution  is  2.5°  longitude  by  4.0°  latitude.  A  90° 
sector  of  the  globe  extending  from  pole  to  pole  is  used  as 
the  computational  domain.  Cyclic  continuity  is  imposed  at 
the  east  and  west  boundaries. 

Convergence  of  the  meridians  toward  the  poles  normally 
requires  a  very  short  time  step  to  ensure  computational 
stability.   Instead  of  shortening  the  time  step,  longitu- 
dinal averaging  is  done  on  the  terms  governing  gravity  wave 
propagation  (Arakawa  and  Lamb,  1977).   The  averaging  is  done 
by  expanding  the  zonal  pressure  gradient  and  zonal  mass  flux 
into  a  Fourier  series  around  each  latitude  circle.   The 
amplitude  of  each  wave  component  is  then  reduced  by  the 
factor  S.   The  gradients  of  the  meteorological  fields  are 
smoothed  instead  of  their  Fourier  expansions.   Application 
of  the  smoothing  operator,  S,  occurs  at  high  latitudes  where 
S  <  1.   The  smoothing  operator  is  defined  as: 


AA_    cos  (j)  ,~  ,, 

Acf>  sin(mAA/2)  {         ' 


(A  list  of  symbol  definitions  is  contained  in  Appendix  B.) 

B.   THE  VERTICAL  GRID 

The  model  equations  are  written  using  the  vertical  sigma 
coordinate: 


16 


p-pt 


a   =  — n —  (2.2) 


where : 


n  =  ps  -  Pt- 

p   is  the  surface  pressure,  and  p.  is  the  pressure  at  the 
top  of  the  model  atmosphere  (p   =  200  mb) .   The  variables 
are  vertically  staggered  in  six  equally  spaced  sigma  layers 
(Figure  2.2).   The  surface  and  terrain  pressure  are  defined 


p  =  200mb  a  =  0 cr=0 

u.v.e 

a 

u.v.e 

— . a 

uv,e 

, . CJ 

— — -  a 

—  U.V0 

a 

p=Ps  CT  =  1  CT=Q 


Figure  2.2.   Vertical  distribution  of  variables  in  a 

six-layer  sigma  coordinate  system.   Solid 
lines  are  a-levels;  dashed  lines  represent 
the  middle  of  the  layer  (Arakawa  and 
Lamb,  1977)  . 


on  the  lower  boundary,  o  =  1.   The  vertical  velocity,  o, 
vanishes  at  the  upper  and  lower  boundaries  and  is  diagnostically 
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derived  on  sigma  levels.   Wind  (u,v)  and  potential  tempera- 
ture (0)  are  computed  at  mid-layer. 

C.   FINITE  DIFFERENCING  SCHEMES 

The  model  uses  two  time  differencing  methods — an  Euler- 
backward  scheme  and  the  centered  time  differencing  (or 
leapfrog)  scheme.   The  model  integration  proceeds  in  a 
series  of  one  Eurler-backward  time  step  followed  by  four 
centered  time  steps.   The  Euler-backward  scheme  eliminates 
the  computational  mode  of  the  centered  scheme  and  selectively 
damps  high  frequency  waves.   The  time  step  for  these  model 
integrations  is  360  s. 


18 


III.   MODEL  INITIALIZATION 

An  analytical  initialization  scheme  developed  by  Hayes 
(19  85)  is  used.   In  this  scheme,  the  initial  mean  wind  is 
specified  analytically  using  the  Bickley  jet  described  by 
Haltiner  and  Williams  (1980) .   The  initial  surface  pressure 
is  derived  from  the  wind  using  the  geostrophic  wind  equation, 
and  the  initial  temperature  is  derived  using  the  thermal 
wind  equation.   In  this  research,  a  stable  planetary  wave 
and  a  synoptic  scale  perturbation  are  superimposed  on  a 
horizontally  and  vertically  varying  mean  current. 

A.   INITIALIZATION  OF  THE  MEAN  FLOW 

A  baroclinic  jet  is  located  at  45°  N  and  has  a  maximum 
velocity  of  approximately  46  m/s  at  200  mb .   The  initial 
mean  surface  flow  is  given  by  the  following  equations: 


UC((H   =   usech2  (y  ((})-(},)  )  (3.1) 

o  ill  \J 


v  ((J>)   =   0  .  (3.2) 

s  T 


It  has  the  same  latitudinal  variation  as  the  jet  with  a 
maximum  speed  of  5  m/s  at  45°  N  decreasing  to  zero  at  the 
poles.   The  v-component  of  the  wind  is  initially  set  to 
zero.   The  mean  upper-level  wind  is  given  by: 
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2  £n(ps/p) 


uu(cj>,p(a,4>)  )   =   um  sech  (y(<$>-<$>q)  )  (£n(p  /p  }  )  (3.3) 


s  m 


vu(<}),p(af<t)))=0  (3.4) 


The  wind  varies  in  the  vertical  as  the  natural  logarithm  of 
pressure. 

The  mean  surface  pressure  is  calculated  by  integrating 
the  geostrophic  wind  equation: 


4> 
PSU)   =   exp[£n  PSU0)  -  (a/R)    /   fug  ( $)  /Tg  (<J>)  d<j>]       (3.5) 

*0 


where  p   (<j>Q  =  45°  N)  =  1013.25  mb  and  T  {<$>)    represents  the 
mean  surface  temperature.   Because  the  latitude  structure  of 
T  (<}>)  is  unknown,  p  (<j))  is  computed  iteratively  until  the 
adjustment  of  T   to  the  NACA  standard  for  that  p   is  less 
than  0.01°  K. 

The  mean  surface  temperature  is  calculated  by  integrating 
the  geostrophic  thermal  wind  equation  which  gives: 


u     a2fl  <t> 

T(p,d>)       =      T(pf(J)0)     -    R  9n  (p/p--)  /       sin    <f>  sech     <  <J>  —  <f>  Q )  c3. 4>       (3.6) 


The  integrals  are  computed  using  the  Simpson  Rule  approxima- 
tion.  T(p,(j)Q)  is  obtained  from  the  NACA  standard  atmosphere 
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T(p,cf>0)   =   T(p0,<j>0)  -  TZ  (3.7 


where : 


To      p  rR/g 
p0 


The  value  of  TQ  ,  288°  K,  is  taken  from  the  NACA  standard 
atmosphere  using  a  standard  surface  pressure  of  1013.25  mb . 
(Constants  are  defined  in  Appendix  B.) 

B.   INITIALIZATION  OF  THE  HAURWITZ  WAVE 

The  initialization  of  the  Haurwitz  wave  is  accomplished 
in  a  subroutine  developed  by  Dr.  T.  Rosmond.    The  analyti- 
cal solution  for  a  Rossby  wave  on  the  surface  of  a  sphere 
is  described  by  Haurwitz  (1940). 

The  Haurwitz  wave  has  a  zonal  wave  number  4  (Figure 
3.1) .   The  maximum  meridional  velocity  is  16  m/s .   The  u  and 
v  components  are  given  by  the  following  equations: 


U'   =   -^  (3.8 

a3(p 


a  COS  tj)d  A 


The  equations  of  motion  are  used  to  derive  the  perturbation 
pressure : 


Personal  communication,  present  affiliation  Naval 
Environmental  Prediction  Research  Facility,  Monterey,  CA, 
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Figure  3.1. 


Sea-level  pressure  field  for  a  stationary 
Haurwitz  wave  (3.10).   Maximum  meridional 
velocity  is  16  m/s .   Contour  interval  is 
20  mb. 
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T-*  i      „»,  .«,    ,n.       ^s   rn+l   .    m   n-m+1    m   ,       ,-,■,«» 

P'      C2  (oo+ft)  cos  Bt+mA)  [ sin  (^p   -  — ? — -rr-r-  p  ,  ,  J       (3. 10) 

n      Ttn    n(n+l)  ^n+1 


where  p  (sin  tj>)  is  the  associated  Legendre  polynomial  and  co 
is  the  mean  flow. 

The  equations  listed  above  are  for  a  Haurwitz  wave 
superimposed  on  a  mean  flow.   The  Bickley  jet  and  the 
Haurwitz  wave  schemes  both  have  mean  flow  contributions,  and 
one  of  the  mean  flows  needs  to  be  removed.   Instead  of  using 
the  formal  analytical  approach,  the  mean  u,  v  and  p  for 
each  J  (latitude  strip)  is  subtracted  from  the  Haurwitz  solu- 
tion to  provide  the  Haurwitz  wave  contribution.   This  method 
allows  the  sea-level  pressures  to  be  computed  at  the  same 
time  as  the  winds  for  the  Haurwitz  wave  and  the  Bickley  jet. 
Otherwise,  the  winds  would  have  to  be  computed  and  then  the 
pressures.   The  error  in  this  approach  should  not  have  a 
significant  effect  on  the  results. 

The  Haurwitz  wave  u,  v  and  p  contributions  are  added  to 
the  Bickley  jet  and  mean  surface  flow  values.   The  surface 
values  are  spread  to  the  upper  levels  using  the  equations 
described  in  the  previous  section. 

C.   THE  PERTURBATION 

The  perturbation  is  a  barotropic  disturbance  consisting 
of  a  longitudinal  sine  wave  which  is  used  for  baroclinic 
instability  studies.   The  maximum  latitudinal  amplitude  in 
the  meridional  direction  occurs  at  45°  N.   The  following 
equation  describes  the  perturbation: 
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■   =   fQ  (A  sin  (mA)sirT2(2  4>)  )  (3.11) 


The  hydrostatic  equation  and  the  geostrophic  wind 
equation  are  used  to  derive: 


P§'        pf 0 ,  2 

p'   =  ~   ^§f     =      "  ^(Asin(mA)sin  (2<|>))  (3.12) 


7?-  ^T  =   -  ^-(A  sin  (mA)2sin(4<j>))  (3.13) 

arn  dcp        a 


(m  Acos  (mA)sin2  (2(f))  )     (3.14 


i 


af  0  cos  <J)  9<})      a  cos  cf> 


where  A  =  a-1  m/s.   These  perturbation  quantities  are  added 
to  the  mean  current  and  integrated  forward  in  time  to 
produce  the  forecast. 
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IV.   ANALYSIS  OF  RESULTS 

The  hypothesis  is  that  a  cold  surge  results  from  the 
interaction  of  a  growing  wave  number  8  cyclone  moving  east- 
ward and  a  quasi-stationary  stable  forced  wave  number  4. 
Hayes  (1985)  showed  that  the  basic  baroclinic  current,  which 
is  given  in  Equations  (3.1)  and  (3.3),  is  unstable  for  wave 
number  8  but  is  stable  for  wave  number  4.   The  stationary 
wave  could  be  partially  forced  by  the  Himalayas.   Three 
experiments  were  run  to  verify  the  hypothesis.   Experiment 
I  shows  the  effect  of  the  stable  wave  number  4  on  the  south- 
ward migration  of  cold  air.   The  effect  of  cyclogenesis  in 
wave  number  8  is  depicted  Experiment  II.   Experiment  III 
is  a  combination  of  Experiments  I  and  II. 

A.   EXPERIMENT  I 

In  this  experiment,  the  model  is  initialized  with  the 
stable  wave  number  4  and  the  baroclinic  mean  flow  described 
in  Chapter  III. 

The  baroclinic  jet  is  confined  to  the  northern  hemisphere 
to  minimize  forcing  from  the  southern  hemisphere.   The 
model  is  integrated  forward  in  time  to  60  h.   Sea-level 
pressure  fields  for  12,  24,  36,  48  and  60  h  are  shown  in 
Figures  4.1  to  4.5.   The  potential  temperature  fields  at 
a  =  11/12  for  the  same  times  are  displayed  in  Figures  4.6  to 
4.10.   The  corresponding  wind  vectors  are  shown  in  Figures 
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Figure  4.1 


Sea-level  pressure  field  at  12  h  for 
Experiment  I.   Contour  interval  is  5  mb 
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Figure  4 . 2 


Sea-level  pressure  field  at  24  h  for 
Experiment  I.   Contour  interval  is  5  mb 
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Figure    4.3. 


Sea-level  pressure  field  at  36  h  for 
Experiment  I.   Contour  interval  is  5  mb 
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Figure  4.4. 


Sea-level  pressure  field  at  48  h  for 
Experiment  I.   Contour  interval  is  5  mb 
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Figure    4.5. 


Sea-level  pressure  field  at  60  h  for 
Experiment  I.   Contour  interval  is  5  mb 
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Figure  4.6 


Potential  temperature  field  at  12  h  for 
Experiment  I.   Contour  interval  is  5°  K 
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Figure  4.7. 


Potential  temperature  field  at  24  h  for 
Experiment  I.   Contour  interval  is  5°  K 
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Figure    4.8. 


Potential  temperature  field  at  36  h  for 
Experiment  I.   Contour  interval  is  5°  K, 
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Figure    4.9 


Potential  temperature  field  at  48  h  for 
Experiment  I.   Contour  interval  is  5°  K 
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Figure  4.10. 


Potential  temperature  field  at  60  h  for 
Experiment  I.   Contour  interval  is  5°  K. 
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4.11  to  4.15.   The  90°  sectors  are  plotted  twice  to  produce 
hemispheric  displays. 

In  the  first  twelve  hours,  the  effects  of  the  baroclinic 
jet  on  the  wave  are  readily  apparent  (Figure  4.1).   The  wave 
has  moved  westward  approximately  5°  compared  to  10°  in  the 
southern  hemisphere.   The  westward  movement  is  caused  by  the 
beta  effect.   A  slight  NE-SW  tilt  in  the  sea-level  pressure 
field  is  apparent.   The  absence  of  a  corresponding  tilt 
in  the  southern  hemisphere  is  consistent  with  the  absence  of 
the  jet.   In  the  northern  hemisphere,  the  tilt  is  in  the 
same  sense  as  the  shear  on  the  south  side  of  the  jet  indi- 
cating that  the  wave  will  lose  energy  to  the  mean  flow.   The 
wave  should  damp. 

The  pressure  gradient  has  weakened  as  the  cyclone  and 
anticyclone  are  reduced  in  amplitude  by  5  mb.   However,  at 
the  same  time,  the  amplitude  of  wave  number  4  in  the  poten- 
tial temperature  field  (Figure  4.6)  has  increased.   The 
305°  K  isotherm  has  moved  south  to  near  10°  N.   The  jet  has 
advected  that  field  eastward  to  a  position  45°  out  of  phase 
with  the  pressure  field.   Northeasterly  winds  of  approxi- 
mately 13  m/s  are  observed  at  15°  N  (Figure  4.11). 

By  24  h,  the  tilt  in  the  sea-level  pressure  trough  has 
increased  to  produce  a  stronger  pressure  gradient  in  the 
northwest  quadrant  and  a  weaker  gradient  in  the  northeast 
quadrant  of  the  anticyclone.   There  has  been  a  corresponding 
change  in  the  wind  field  (Figure  4.12) .   A  large  area  of 
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Wind  vectors  at  12  h  for  Experiment 
I.   8  mm  =  20  m/s . 
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Figure  4.12 


Wind  vectors  at  24  h  for  Experi- 
ment I.   7  mm  =  20  m/s. 
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Figure    4.13. 


Wind  vectors  at  36  h  for  Experi- 
ment I.   7  mm  =  20  m/s. 
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Figure    4.14. 


Wind   vectors    at    48    h    for    Experi- 
ment   I.       8    mm   =    20    m/s . 
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Figure  4.15. 


Wind  vectors  at  60  h  for  Experi- 
ment I.   7  mm  =  20  m/s . 
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northeasterly  winds  is  present  in  the  southeast  quadrant. 
Pressure  centers  have  continued  to  weaken.   A  well-defined 
trough  has  formed  along  the  equator  with  wind  speeds  as 
high  as  13  m/s.   The  tilt  and  amplitude  of  the  potential 
temperature  wave  have  increased.   Strong  southwesterlies 
ahead  of  the  trough  are  advecting  warm  air  northward  from 
20°  N  to  45°  N. 

At  36  h,  the  anticyclone  has  intensified  by  about  5  mb. 
Ridging  dominates  the  pattern  particularly  south  of  30°  N. 
The  trough  has  filled  as  the  flow  has  become  more  zonal. 
The  cyclone  has  intensified  slightly.   Weak  cyclonic  flow  is 
evident  at  the  base  of  the  trough  (Figure  4.13)  .   Amplifi- 
cation of  the  potential  temperature  field  has  continued  with 
strong  warm  air  advection  ahead  of  the  trough.   Cold  air 
advection  behind  the  trough  increased  south  of  20°  N.   With 


cold  air  moving  south  and  warm  air  moving  north,  V'0'  is 
positive,  which  indicates  a  baroclinic  process.   The  area 
and  speed  of  northeasterly  winds  on  the  eastern  side  of  the 
anticyclone  have  increased. 

At  48  h,  the  ridging  south  of  30°  N  has  expanded  east- 
ward while  the  anticyclone  weakened.   North  of  45°  N  the 
cyclone  has  intensified  slightly  and  enlarged.   The  pressure 
gradient  is  oriented  north-south.   Although  not  evident  on 
the  pressure  field,  a  weak  cyclone  has  formed  at  the  base  of 
the  trough  (Figures  4.4  and  4.14).   The  wave  structure  in  the 
potential  temperature  has  amplified  further.   The  305°  K 
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isotherm  has  moved  to  within  3°  of  the  equator.   North  of 
30°  N,  the  field  is  180°  out  of  phase  with  that  south  of 
30°  N.   Winds  are  northeasterly  on  the  upwind  side  of  the 
thermal  trough. 

In  the  next  twelve  hours,  rapid  development  has  occurred 
in  the  cyclone.   The  central  pressure  has  dropped  almost 
10  mb.   A  secondary  anticyclone  formed  15°  northwest  of  the 
primary  high  pressure  center.   The  wave  in  the  potential 
temperature  field  has  continued  to  amplify.   The  305°  K 
isotherm  has  reached  1°  S.   The  field  does  not  reflect  the 
surface  pressure  pattern  north  of  30°  N.   Northeasterly 
winds  of  greater  than  15  m/s  on  the  eastern  side  of  the 
primary  anticyclone  are  observed  in  the  equatorial  region 
(Figure  4.15).   The  secondary  anticyclone  is  not  adequately 
reflected  in  the  wind  fields  indicating  that  this  feature 
is  not  in  geostrophic  balance. 

This  experiment  demonstrates  the  effects  of  a  baroclinic 
jet  on  the  stable  wave  number  4.   The  jet  provided  enough 
forcing  to  produce  cold  temperatures  and  strong  northerly 
winds  in  the  tropics.   The  305°  K  isotherm  moved  southward 
from  near  10°  N  to  south  of  the  equator.   At  60  h,  the 
secondary  pressure  centers  indicate  growth  of  wave  number  8 
which  is  unstable.   This  wave  is  induced  by  the  nonlinear 
interaction  of  the  long  wave  with  itself. 
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B.   EXPERIMENT  II 

In  this  experiment,  the  model  was  initialized  with  the 
baroclinically  unstable  mean  flow  and  the  wave  number  8 
disturbance  described  in  Chapter  III.   The  model  was 
again  integrated  forward  in  time  to  60  h.   Sea-level 
pressure  fields  for  12,  24,  36  and  48  h  are  shown  in  Figures 
4.16  through  4.19.   The  potential  temperature  fields  for 
12,  24,  36  and  48  h  are  depicted  in  Figures  4.20  through  4.23 

The  results  of  this  experiment  are  used  to  determine 
when  to  add  the  stable  wave  number  4  flow  which  will  advect 
the  synoptic  scale  disturbance  southward  in  Experiment  III. 
One  requirement  for  this  simulation  is  to  have  a  sea-level 
pressure  difference  of  at  least  25  mb  between  high  and  low 
centers  to  ensure  an  adequate  northwesterly  current  for 
surge  initiation  (Baker,  1983).   Another  requirement  is  the 
appearance  of  a  high  pressure  center  representing  the 
Siberian  anticyclone  to  the  northwest  of  the  low  pressure 
center. 

As  shown  in  Figure  4.16,  the  disturbance  is  primarily  a 
feature  of  the  northern  hemisphere  mid-  and  upper-latitudes 
with  a  very  flat  pressure  gradient  in  the  tropical  region. 
In  the  first  twelve  hours  a  low  has  formed  in  the  trough 
approximately  at  60°  N.   The  trough  has  a  slight  NE-SW 
tilt  and  has  moved  eastward  approximately  10°  at  45°  N. 
A  sea-level  pressure  difference  of  greater  than  25  mb  exists 
between  high  and  low  centers.   The  potential  temperature  and 
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Figure  4.16 


Sea-level  pressure  field  at  12  h  for 
Experiment  II.   Contour  interval  is  5  mb 
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Figure  4 . 17 


Sea-level  pressure  field  at  24  h  for 
Experiment  II.   Contour  interval  is  5  mb 
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Figure  4.18. 


Sea-level  pressure  field  at  36  h  for 
Experiment  II.   Contour  interval  is  5  mb 


47 


30    N 


^y/KM^-^Mi'   \\r'  '^-v?: ■-■,'•••■    ^~-v:  "  -   ,  .1- ?=-- -;^< &• x 

00W§§  c  jm  C  f9 


60 


/  v  f  VS- 


\    X-^A 


Figure    4.19. 


Sea-level  pressure  field  at  48  h  for 
Experiment  II.   Contour  interval  is  5  mb 
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Figure    4.20 


Potential  temperature  field  at  12  h  for 
Experiment  II.   Contour  interval  is  5°  K 
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Figure    4.21. 


Potential  temperature  field  at  24  h  for 
Experiment  II.   Contour  interval  is  5°  K 
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Figure    4.22. 


Potential  temperature  field  at  36  h  for 
Experiment  II.   Contour  interval  is  5°  K 
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Figure    4.23. 


Potential  temperature  field  at  48  h  for 
Experiment  II.   Contour  interval  is  5°  K, 
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wind  fields  indicate  a  baroclinic  process  with  warm  air 
being  advected  northward  and  cold  air  being  advected 
southward. 

At  24  h,  the  low  centers  have  moved  5°  southeast  of 
their  previous  positions  and  have  intensified.   A  sharp 
surface-pressure  ridge  extending  north  to  80°  N  separates 
the  disturbances .   Eastward  movement  of  the  trough  has 
continued  at  the  previous  rate.   Development  of  the  primary 
disturbance  in  the  wind  and  potential  temperature  fields 
parallels  that  of  the  pressure  field.   Some  short  wave  noise 
is  apparent  in  the  wind  field.   Noise  is  to  be  expected 
because  no  surface  frictional  dissipation,  diffusion  or 
smoothing  techniques  are  included  in  the  model.   These 
techniques  tend  to  dissipate  surge  energy,  and  the  objective 
of  this  study  is  the  simulation  of.  intense  monsoon  surges. 

From  24  h  to  36  h,  rapid  deepening  has  occurred  in  the 
cyclone.   The  central  pressure  of  the  low  has  dropped  at 
least  15  mb.   The  ridge  between  low  centers  has  weakened, 
but  the  sea-level  pressure  difference  between  high  and  low 
centers  has  increased  to  more  than  50  mb.   A  broad  band  of 
low  pressure  appears  to  dominate  the  zone  between  50°  N  and 
60°  N.   Surface  warm  and  cold  fronts  are  clearly  evident 
in  the  temperature  field  at  36  h  (Figure  4.22).   The  potential 
temperature  shows  a  poleward  flux  of  warm  air  and  a  migration 
of  cold  air  south  of  30°  N. 
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The  rapid  deepening  of  the  surface  pressure  low  con- 
tinued between  36  h  and  48  h.   The  central  pressure  dropped 
approximately  10  mb.   High  pressure  centers  formed  near 
40°  N,  and  the  sea-level  pressure  difference  between  the 
high  and  low  pressure  centers  is  greater  than  70  mb.   The 
northern  portion  of  the  strong  baroclinic  zone  is  advected 
southeastward  about  7°  latitude.   The  potential  temperature 
field  shows  some  additional  migration  of  cold  air  south  of 
30°  N  (Figure  4.23)  compared  to  12  h  earlier.   At  60  h,  the 
300°  K  isotherm  follows  the  30°  N  latitude  circle,  and  the 
305°  K  isotherm  is  a  little  south  of  20°  N. 

Analysis  of  these  fields  reveals  several  key  features. 
A  sea-level  pressure  difference  of  greater  than  25  mb 
between  high  and  low  centers  is  present  from  12  h  onward. 
The  sea-level  pressure  difference  actually  increases 
throughout  the  run.   A  high  pressure  center  northwest  of  the 
low  is  best  depicted  at  24  h;  it  weakens  afterwards.   A 
strong  Siberian  high  is  not  observed  probably  due  to  the 
presence  of  a  strong  jet  and  lack  of  topographic  and  radia- 
tive forcing.   The  period  of  most  rapid  cyclonic  develop- 
ment is  between  24  h  and  36  h.   The  cold  air  and  strong 
northerly  winds  never  penetrate  south  of  30°  N. 

Based  on  these  results,  wave  number  4  should  be  added  to 
the  growing  cyclone  between  24  h  and  36  h.   The  most  rapid 
cyclogenesis  is  occurring  at  this  point.   Adding  the  larger 
scale  wave  will  further  advect  the  cold  air  south  to  help 
initiate  the  surge. 
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C.   EXPERIMENT  III 

The  initialization  for  this  experiment  is  the  same  analy- 
tic scheme  as  Experiment  II.   However,  after  36  h  of  model 
time,  the  wave  number  4  fields  are  added.   (Note:   The 
zonal  mean  has  been  removed  from  the  wave  number  4  fields.) 
Sea-level  pressure,  potential  temperature,  and  wind  fields 
are  depicted  in  Figures  4.24  through  4.35. 

The  development  of  the  disturbance  through  36  h  is 
described  in  Experiment  II.   The  effect  of  wave  number  4  on 
the  disturbance  is  readily  apparent  when  comparing  the  36  h 
fields  of  Experiments  II  and  III.   A  high  pressure  center 
of  at  least  1030  mb  forms  in  the  ridge.   The  primary  cyclone 
deepens  another  35  mb  with  the  addition  of  wave  number  4. 
This  rapid  deepening  produces  a  pressure-  surge . 

The  second  stage  of  the  monsoon  surge  is  characterized 
by  a  sharp  drop  in  surface  dewpoint  (Baker,  1983).   This 
model  does  not  include  moisture  or  its  effects,  but  observa- 
tional studies  (Chang  et  al . ,  1979  and  Chu,  1978)  show  that 
this  stage  can  also  be  associated  with  a  drop  in  surface 
temperature.   At  36  h,  a  trough  of  cold  air  extends  south- 
ward to  approximately  17°  N  at  the  base  of  the  high  pressure, 
however,  the  potential  temperature  gradients  are  very  weak 
south  of  30°  N.   The  winds  are  primarily  northerly  at 
15  m/s  on  the  eastern  side  of  the  anticyclone. 

From  36  h  to  48  h  the  weaker  low  moved  east  toward  the 
stronger  low.   The  long  wave  ridge  weakened  and  tilted  more 
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Figure  4.24 


Sea-level  pressure  field  at  36  h  for 
Experiment  III.  Contour  interval  is 
10  mb. 
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Figure    4.25. 


Sea-level  pressure  field  at  48  h  for 
Experiment  III.  Contour  interval  is 
10  mb. 
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Figure    4.26. 


Sea-level  pressure  field  at  54  h  for 
Experiment  III.  Contour  interval  is 
10  mb. 
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Figure    4.27. 


Sea-level  pressure  field  at  66  h  for 
Experiment  III.  Contour  interval  is 
10  mb. 
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Figure    4".  2 


Potential  temperature  field  at  36  h  for 
Experiment  III.   Contour  interval  is 
5°  K. 
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Figure  4.29. 


Potential  temperature  field  at  48  h  for 
Experiment  III.   Contour  interval  is 
5°  K. 
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Figure  4-30. 


Potential  temperature  field  at  54  h  for 
Experiment  III.   Contour  interval  is 
5°  K. 
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Figure  4.31. 


Potential  temperature  field  at  66  h  for 
Experiment  III.   Contour  interval  is 
5°  K. 
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Figure    4 . 32 . 


Wind   vectors    at    36    h    for   Experiment 
III.       4    mm   =    20   m/s. 
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Figure    4.33 


Wind  vectors  at  48  h  for  Experiment 
III.   8  mm  =  50  m/s . 
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Figure  4.34. 


Wind  vectors  at  54  h  for  Experiment 
III.   8  nun  =  50  m/s. 
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Figure    4.35. 


Wind  vectors  at  66  h  for  Experiment 
III.   8  mm  =  50  m/s . 
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NE-SW.   The  small  cyclone  deepened  at  least  10  mb .   The 
large  cyclone  deepened  approximately  10  mb  and  moved 
slightly  eastward  with  the  mean  flow.   The  tilting  and 
weakening  of  the  ridge  caused  three  anticyclonic  circulations 
to  develop  within  the  ridge  particularly  toward  the  north- 
west.  At  4  8  h,  winds  are  northeasterly  from  40°  N  to  about 
10°  N  with  a  speed  of  at  least  18  m/s. 

The  48  h  potential  temperature  field  shows  a  broad 
tongue  of  cold  air  extending  southwest  near  20°  N  as  if 
spreading  over  the  South  China  Sea.   The  furthest  extent  is 
approximately  15°  N.   The  thermal  gradient  remains  weak  south 
of  30°  N  and  especially  south  of  20°  N.   A  strong  baroclinic 
zone  is  depicted  behind  the  cold  front. 

Six  hours  later  the  secondary  low  pressure  center  is 
located  5°  latitude  further  southeastward.   The  primary  low 
center  remained  stationary.   The  ridge  elongated  and  moved 
westward.   This  movement  resulted  in  a  band  of  northeasterly 
winds  (of  approximately  20  m/s)  extending  southwestward 
from  40°  N  to  at  least  10°  N  (Figure  4.34)  as  if  along  the 
eastern  coast  of  China.   These  northeasterly  winds  advect 
colder  temperatures  as  far  south  as  15°  N  (Figure  4.30) . 
The  strong  baroclinic  zone  moved  approximately  2.5° 
southeastward. 

From  54  h  to  66  h  the  low  pressure  center  moves  east- 
ward, and  the  high  pressure  center  moves  westward.   The 
pressure  gradient  south  of  30°  N  decreases.   The  NE-SW 
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band  of  northeasterly  winds  narrowed  but  continues  to 
advect  cold  temperatures  as  far  south  as  approximately 
10°  N.   The  strong  baroclinic  zone  continued  moving  south- 
eastward but  at  a  slower  rate. 

Several  key  observations  are  made  based  on  this  analy- 
sis.  Northeasterly  winds  move  rapidly  south  to  the  tropical 
regions  in  12  h.   These  winds  brought  cold  mid-latitude 
air  to  the  tropics  in  30  h.   Compared  to  Experiment  I,  the 
northeasterly  winds  are  stronger  but  do  not  go  quite  as 
far  south  into  the  equatorial  region.   The  sea-level  pressure 
pattern  progresses  eastward  very  little  because  additional 
forcing  is  necessary  to  overcome  the  beta  effect  on  wave 
number  4.   The  pressure  ridge  and  the  cold  air  tongue  are 
better  defined  than  the  wind  surge.   Baroclinic  processes 
are  present  in  the  results  of  this  experiment  south  of 
30°  N. 


69 


V.   CONCLUSIONS 

The  interaction  of  a  stable  planetary  wave  with  a 
growing  cyclonic  disturbance  in  a  mid-latitude  baroclinic 
flow  was  examined.   The  purpose  of  this  study  was  to  deter- 
mine the  contribution  of  a  quasi-stationary  stable  wave 
number  4  to  the  initiation  of  a  monsoon  surge  as  a  wave 
number  8  cyclone  develops. 

Several  circulation  features  have  been  identified  as 
characteristic  of  monsoon  surges.   They  include:   (1)  the 
strengthening  of  north-northeasterly  winds  in  the  South 
China  Sea,  (2)  a  sharp  drop  in  surface  temperature  at  obser- 
vation stations  in  China  and  Southeast  Asia,  (3)  the  shallow 
structure  of  the  surge,  (4)  the  strong  baroclinic  zone  and 
(5)  northerly  winds  in  the  tropics.   In  varying  degrees, 
each  of  these  characteristics  were  observed  in  one  or  more 
of  the  experiments . 

In  Experiment  I,  the  jet  provided  enough  additional 
forcing  to  the  wave  number  4  to  produce  cold  temperatures 
and  strong  northerly  winds  in  the  tropics.   The  305°  K 
isotherm  moved  from  near  10°  N  to  south  of  the  equator  in 
48  h.   The  winds  were  initially  northerly  with  a  small  area 
of  northeasterly  winds.   Throughout  the  experiment,  the  winds 
on  the  east  side  of  the  anticyclone  became  more  northeasterly 
and  increased  in  speed  in  the  equatorial  region.   A  strong 
baroclinic  zone  was  observed  on  the  downwind  side  of  the  trough 
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Experiment  II  showed  the  evolution  of  an  unstable  baro- 
clinic  wave  which  was  limited  to  the  mid-latitudes.   The 
300°  K  isotherm  moved  south  approximately  5°  latitude  in 
48  h,  as  opposed  to  10°  and  south  of  the  equator  in  Experi- 
ment I.   A  strong  baroclinic  zone,  which  is  characteristic 
of  a  frontal  zone,  was  evident  in  the  low-level  potential 
temperature  fields  at  the  mid-latitudes,  but  it  did  not 
propagate  southward  into  the  tropics.   In  this  experiment 
there  were  no  significant  northerly  winds  south  of  45°  N 
and  no  cold  temperatures  south  of  30°  N. 

Several  characteristics  of  monsoon  surges  were  observed 
in  Experiment  III  in  which  the  wave  number  4  and  wave  number 
8  were  combined.   North-northeasterly  winds  reached  the 
tropical  regions  12  h  after  the  addition  of  the  long  wave. 
Cold  air  moved  southward  to  approximately  10°  N  within  30  h 
after  the  addition  of  the  long  wave.   A  strong  baroclinic 
zone  was  evident  in  the  potential  temperature  fields.   It 
was  advected  southward  to  the  tropics.   Northerly  winds  were 
observed  in  the  tropics  at  least  as  far  south  as  10°  N. 
Experiment  III  roughly  approximated  events  occurring  during 
a  cold  surge. 

These  results  suggest  that  the  combination  of  quasi- 
stationary  stable  long  wave  and  a  developing  mid-latitude 
cyclone  initiate  a  monsoon  surge.   Additional  forcing  appears 
to  be  required  to  keep  wave  number  4  stationary.   Wave  number 
4  also  should  be  modified  so  that  it  does  not  extend  to  the 
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equator  initially.   Moving  the  jet  and  the  maximum  lati- 
tudinal amplitude  of  the  perturbation  to  35°  N  would  increase 
the  baroclinicity  south  of  30°  N  as  well  as  more  closely 
model  the  observed  position  of  the  jet  during  a  surge  event. 
This  study  found  that  a  quasi-stationary  stable  wave  number 
4  and  an  unstable  wave  number  8  could  initiate  a  monsoon 
surge . 
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APPENDIX  A 
MODEL  PRIMITIVE  EQUATIONS 


Equation  of  State: 


pa   =   RT 


Hydrostatic  Equation: 


-  it  a 


Horizontal  Momentum  Equation: 


n(g)  +  f  kn  v  +  v0(no)  -  v  n^-^-  =  o 


Continuity  Equation: 


M+v0-,nv,  +i|f>  -  o 


1st  Law  of  Thermodynamics: 


3(nC  £n9)         ~  9(naC   £n6) 

+  v-(nvc  in  e)  +  ~-£ =  o 


3a  p  3a 
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APPENDIX  B 
LIST  OF  SYMBOLS 

A  constant 

a  radius  of  the  earth 

C  constant 

C  specific  heat  at  constant  pressure 

f  Coriolis  parameter 

ffl  value  of  Coriolis  parameter  at  45°  N 

g  gravitational  acceleration 

m  wave  number 

p  pressure  (p  =  p  +  p'  +  P1) 

p0  surface  pressure  at  45°  N,  1013.25  mb 

p  maximum  model  pressure  height,  12  5  mb 

p  surface  pressure 

p  pressure  at  tropopause 

R  gas  constant  for  air 

S  smoothing  operator 

T  temperature 

TQ  surface  temperature  at  45°  N,  288°  K 

T  surface  temperature 

u  zonal  velocity  component  (u  =  u  +  u1  +  U1) 

u  maximum  velocity,  surface:   5  m/s,  upper:   6  5  m/s 

u  surface  velocity  component 

u  upper  level  velocity  component 
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v  meridional  velocity  component  (v  =  v  +  v1  +  V) 

Z  terrain  height 

a  specific  volume 

3  frequency 

y  halfwidth  of  jet 

r  dry  adiabatic  lapse  rate,  6.5°  K/km 

6  potential  temperature 

A  longitude 

tt  3.14159 

II  normalized  pressure 

a  vertical  coordinate 

a  vertical  velocity  in  sigma  coordinate  system 

tj)  latitude 

(})q  latitude  of  jet 

$  geopotential ,  also  current  terrain  height 

0,  angular  velocity  of  the  earth 

oj  angular  velocity  of  air  relative  to  the  earth 
(mean  flow) 

\b  stream  function 
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